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ABSTRACT 

The results of three quasi-objective techniques for 
enhancing the accuracy of 24-hour cloud cover forecasts 
in sparse data areas, using satellite data and numerical- 
ly-produced analyses and prognoses, are presented. All 
three of the techniques involve the modification of the 
analyzed 500-mb relative vorticity, as based on the sat- 
ellite-observed cloud patterns and the associated 
changes to the prognostic fields of 500-mb relative vor- 
ticity, The final 24-hour cloud-cover forecast is based 
on the manually-modified relative vorticity pattern and 
three different advecting currents. These advecting 
currents are based on the prognostic 500-mb height field, 
modified and unmodified, and the prognostic 500-mb space 
mean field. The study was conducted over a large area 
of the North Pacific Ocean for the period 9 March to 
13 April 1972, 
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INTRODUCTION 



I . 

The objective of this study was to investigate the 
possibility of increasing the accuracy of estisting cloud 
cover forecasting techniques in sparse-data areas by us- 
ing satellite data in conjunction with numerical analyses 
and prognoses. In particular, an improvement was sought 
on Bittner’s [1] quasi-objective method of using the ad- 
vection of relative vorticity on the 500-mb surface to 
determine the primary areas of cloud cover. 

Bittner’s quasi-objective method uses current satel- 
lite-observed cloud patterns to modify the analyzed pat- 
terns of relative vorticity. The difference between the 
Initial and modified relative vorticity patterns is then 
subjectively applied to the 24- or 36-hour prognostic 
patterns of relative vorticity. The primary areas of 
cloud cover are then related to the areas of positive 
relative vorticity advect ion (PVA) as determined by the 
unmodified 24- or 36-hour prosnostic 500-mb height 
fields (Method 1). 

Theoretically, an improvement in the cloud cover fore- 
casts could be made if the prognostic 500-mb height 
fields were modified by the amounts associated with the 
changes in the prognostic relative vorticity fields prior 
to determining the areas of positive relative vorticity 
advection (Method 2) . 
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Also, the theoretical development of Fjortoft [2] in 
manually producing a barotropic 24-hour forecast of the 
500-mb height values indicates that using prognostic 
space-mean fields to determine the areas of positive 
relative vorticity advection might show an improvement 
in the verification of 24-hour cloud cover forecasts 
(Method 3) . 

Therefore, this study was specifically undertaken to 
determine which of the above three methods are the most 
accurate for forecasting 24-hour synoptic-scaled cloud 
patterns. All three techniques were used simultaneously 
to produce cloud-cover forecasts over the North Pacific 
Ocean from 20°N to 50°N and from 170^E to 125‘^W. The 
forecast period extended from 9 March to 13 April 1972. 

Modifying the 500-mb height prognoses in addition to 
the relative vorticity prognoses took slightly longer 
than the other two methods, but the resultant improvement 
in the forecasts seemed to indicate that the small amount 
of additional work was well worth the effort involved. 
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II. BACKGROUND 



The oceanic areas of the world are relatively void 
of conventional meteorological observations, especially 
above the earth’s surface. Over the middle-latitude 
areas the synoptic-scaled cloud patterns attain their 
best development. Fortunately, these cloud patterns, 
which are so revealing of synoptic systems to trained 
meteorologists, are easily observed with satellites. 
Consequently, for the past several years, meteorologists 
have striven to incorporate the vast quantity of synoptic- 
scaled information obtained from satellites into useful 
models for improved forecasts. Further, there has been 
considerable effort to incorporate this satellite-derived 
data into numerical analysis schemes, as Nagle and Hayden[7], 
Nagle and Clark [6], Mantei and Workman [5], Zeigler (8] 
and Glaes [3]. Recently Bittner [1] has illustrated a 
method of using these satellite-observed synoptic-scaled 
cloud patterns to manually improve the numerical analysis 
of the 500-mb relative vorticity field. He then carried 
this improvement: forward 24 to 36 hours by manually adjust- 
ing the numerically-produced prognostic 500-mb relative 
vorticity field. Finally, he converted this modified prog- 
nostic field into cloud and precipitation forecasts. 

In an equivalent barotropic numerical model of the 
atmosphere it may be easily shown that the sign of the 
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advection of relative vorticity at the level of nondiver- 
gence (LND) is the same as the sign of the vertical motion 
through this level. In fact, in this simple model of the 
atmosphere the sign of the vertical motion is the same 
throughout a given column of air, if terrain effects are 
neglected. Therefore, positive advection of relative vor- 
ticity at the LND produces positive (upward) vertical 
motion throughout the column of air. In the real atmosphere 
the sign of the advection of relative vorticity, or geo- 
strophic relative vorticity, on the 500-mb pressure surface 
is highly correlated with the sign of the synoptic-scaled 
vertical motion through this surface. That is, if the 
advection of geostrophic relative vorticity on the 500-mb 
surface is negative then the synoptic-scaled vertical motion 
in the vicinity of the 500-mb level is also negative (i.e. 
subsidence) almost all of the time. 

If hygroscopic nuclei and moisture are present within 
the atmosphere then the formation and dissipation of 
clouds are largely determined by the sign of the vertical 
motion. That is, with the above conditions, positive 
vertical motion on a synoptic scale will produce synoptic- 
scaled cloud patterns when the time scale is sufficient 
to cover the phase lag between initial lifting and subse- 
quent condensation. With a given moisture distribution, 
the stronger the advection of vorticity the shorter the 
time lag between the vertical ascent and its tell-tale 
signature, the synoptic-scaled cloud patterns. 
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For a given analysis or prognosis of the height field 
of a constant pressure surface the relative vorticity may 
be obtained diagnostically from an equation which relates 
the mass field to the wind field. The exact form of this 
diagnostic equation depends upon the assumptions made 
concerning the relationship between the wind field and 
the mass field. Fjortoft [2] demonstrated that the height 
field of a constant pressure surface may be smoothed by 
using a space differential operator called the Laplacian 
or del-squared operator. This smoothed height field is 
more commonly known in meteorology as the space mean. The 
space mean represents the initial height field after the 
shorter wave lengths have been selectively removed. The 
exact short wave length which is completely removed from 
this field depends upon the grid length selected for the 
finite difference equation used in the approximation of 
the analytical equation. Fjortoft proved that, except for 
a latitude-dependent weighting function, the geostrophic 
relative vorticity field is obtained when the initial 
height field is subtracted from this space-mean height 
field. 

The height field of a constant pressure surface which 
is obtained through numerical analysis or prognosis may 
be subjected to the scale separation technique of Roll [4]. 
In effect the initial height field is repeatedly smoothed 
with a Laplacian operator until a residual field called 
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the SR field remains. This SR field may be roughly con- 
strued to be equivalent to the space-mean field described 
above. When this SR field is subtracted from the initial 
height only a disturbance field, called an SD field, re- 
mains. The SD field roughly represents the combined effect 
of all the numerically representable short waves with wave 
lengths shorter than wave number six. Likewise the SR pat- 
tern approximately depicts wave numbers zero (polar vortex) 
through six. Note that the field which approximates the 
geostrophic relative vorticity is the space-mean minus the 
initial height field while the SD field is minus the SR 
(space mean) plus the initial height field. That is, when 
the pattern of the SD analysis is used to approximate the 
pattern of the geostrophic relative vorticity, the sign of 
the SD centers must be reversed. See Appendix A. 

Consequently, the above concepts led to the development 
of quasi-objective rules for relating synoptic-scale clouds 
to the 500-mb SD field, Bittner [1] outlined these and for 
convenience they are listed in Appendix B. Bittner’s 
method consists of using the most recent satellite-observed 
cloud patterns to modify the current 500-mb SD analysis. 

The resulting difference field between the initial and modi- 
fied SD patterns is applied to the prognostic SD field to 
produce modified forecast SD patterns. 

The particular step of modifying the SD patterns of the 
prognostic field based on the modifications applied to the 
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current SD patterns appears to be the most subjective of 
all the procedural steps of Bittner’s method. If there 
is no energy transformation occurring between the short 
waves represented by the SD field and the long waves por- 
trayed by the SR field during the period of the prognosis, 
then the modifications of the prognostic SD patterns are 
rather straight forward. However, for a given cyclonic 
or anticyclonic system within an initial SD field, the 
greater the forecast change of development of a short-wave 
feature the less reliable the modification which is applied 
to the SD system of the prognostic field. That is, the 
more an SD center changes during a prognostic period, the 
more difficult it becomes to propagate an initial modifi- 
cation of the SD center to its prognostic counterpart. 

The modified SD prognosis is then used in conjunction 
with the unchanged 500-mb height prognosis to derive a 
subjective cloud and precipitation forecast valid at the 
time of the numerical height forecast. Bittner makes the 
following assumptions relative to his synoptic-scaled 
cloud pattern forecasts: 

1) Large-scale vertical motion fields are highly cor- 
related with the characteristic distribution of extratropi- 
cal, synoptic-scale cloud patterns as seen by satellites. 

2) In turn, the vertical motion patterns are closely 
related to the advection of relative vorticity at 500 mb. 

3) An SD field is equivalent to a relative vorticity 
field. See Appendix A. 
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4) Any errors in the SD field based on the initial 
500-mb analysis are propagated during the numerical prog- 
nosis to produce similar errors in the SD field based on 
the prognostic 500-mb heights. Hence the error field is 
assumed to be conserved during the prognosis. 
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Ill, DATA SOURCES AND SELECTION OF STUDY AREA AND PERIOD 



An area of the North Pacific Ocean extending from 20°N 
to 50°N and from 170^E to 125°W was selected for investiga- 
tion for several reasons. 

This area was selected primarily because it is both 
analyzed numerically by FNV7C and has almost continuous day- 
light picture coverage by the ATS-1 satellite. ATS-1 is 
the designator for the geostationary Applications Technology 
Satellite positioned near the equator of the Pacific Ocean 
since November 1967. Its spin-scan cloud camera takes time- 
lapse photographs approximately every 25 minutes. The re- 
sulting pictures cover a region from about 50®S to 50°N over 
a span of 100 degrees of longitude centered at 150°W, with 
a resolution of two miles at the center. 

The subject ATS-1 photographs were generally available 
locally from FNWC on a real-time basis Monday through Friday. 
Since the Saturday photograph was needed for verification 
of Friday’s forecast, the study period was usually limited 
to four days per week, namely Monday through Thursday. Some 
of these photographs and others that were missing on occa- 
sion were supplied by the Naval Air Systems Command’s Pro- 
ject FAMOS.^ The FNWC computer products were also available 
locally on a real-time basis. 

^Presently part of the Environmental Prediction Research 
Facility, Monterey, California. 
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The period selected (March and April 1972) not only af- 
forded a real-time basis for investigation but it is also 
the time of year when large, well-defined systems charac- 
terize the North Pacific Ocean, These systems readily lend 
themselves to the modification techniques evaluated during 
this study. 
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IV, PROCEDURES 



The research objectives of the procedures described 
below were to extend Bittner^s work [1] and to outline an 
approach which would be practical for Fleet use. At the 
time when the FNWC 0000 GMT numerical analyses are received 
by the Weather Centrals and Facilities, the prognostic 
fields are not yet available. Therefore, to expedite 24- 
hour cloud forecasts from 0000 GMT it would be best to use 
the 36-hour prognostic height fields, as based on the pre- 
vious 12-hour analysis, and the current 0000 GMT analysis. 
Bittner indicated his approach is good for 24- and 36-hour 
forecasts. Since he effectively assumes that there is no 
energy transformation occurring betv/een the SD and SR 
fields for 24 to 36 hours, then a change in an SD pattern 
may be propagated forward from any time step within this 
forecast period. That is, a change in the SD pattern at 
plus 12 hours fromthe time of the forecast may be carried 
forward to plus 36 hours. Therefore, in these procedures, 
changes noted in the 0000 GMT analysis of the SD patterns 
were incorporated in the prognostic SD patterns which 
verified 24 hours later, but which were made from the 
previous 1200 GMT analysis. 

Another point which may be mentioned in this regard is 
that within a sparse upper-air data area such as the North 
Pacific Ocean, the difference between the 12-hour prognosis 
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of the 500-mb level and the verifying 500-mb analysis may 
be quite small# 

The procedures in this study deal almost exclusively 
with FNWC ’ s scale separation products reproduced on a 
1:60,000,000 polar stereographic projection. This required 
one scale change — that of manually transferring the cloud 
features of the ATS-1 photograph to this scale. FNWC is 
currently testing a program that will convert their pro- 
ducts to the scale and projection of the ATS-1 photograph. 
When this is perfected, all future studies can be carried 
out On the ATS-1 scale to eliminate any errors in the 
manual transfer process. 

In general, this study consisted of determining a modi- 
fied SD field and projecting these modifications on the 
prognostic SD field. Then by using three advecting fields 
(FNWC’s Z and SR PROGS, and a hand-modified FNWC Z PROG, 
called Z’PROG), three cloud cover forecasts were determined 
for verification. The detailed procedures are as follows: 

1) Each day the ATS-1 photograph closest in time to 
0000 GMT was selected for use. The major cloud features 
were traced onto a matching grid (same scale and subpoint 
as the photograph), carefully delineating as covered those 
areas with greater than 50% cloud cover and as open those 
areas with less than 50% coverage. 

2) The cloud patterns were then transferred by hand to 
a polar-stereographic projection of 1:60,000,000 scale. 
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This was a difficult and tedious procedure at first but 
after some practice, the time was decreased to less than 
15 minutes per map, 

3) FNWC’s SD and SR analyses were put in register 
with the transferred cloud patterns. In areas of major 
cloud features such as fronts, cyclones, etc,, the SD 
patterns and isolines were modified in accordance with 

the rules specified by Bittner [1]. These rules are listed 
in Appendix B for convenience to the reader. Further, areas 
of greater than 50% cloud cover were assumed to be associa- 
ted with positive vertical motion as indicated by positive 
vorticity advection using analyzed SR winds as the advecting 
field. The SD lines were redrawn to indicate the required 
advection in order to support the observed cloud patterns. 
Since the zero lines are not produced on the SD charts, 
these lines were drawn by hand. The result of the above w*:s 
a modified SD chart, hereafter referred to as the SD’ 
analysis. The SR fields were not modified. 

4) Subjectively, the same adjustments that were made 
to the SD analysis were then applied to the 36-hour SD 
prognosis, hereafter designated as the SD’ prognosis. Since 
the location of zero-valued SD isolines are somewhat arbi- 
trary, care must be used to surround the same centers with 
the same zero line on both the analysis and the prognosis, 

5) The SD’ prognosis was then graphically added to 
the SR prognosis to produce a modified 500-mb height prog- 
nosis, hereafter indicated as a Z’ prognosis. 
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6) The Z’ prognosis was put in register with the SD’ 
prognosis and used as the advecting field. Clouds were 
forecast in all areas of negative SD advection (positive 
relative vorticity advection). Clear skies (less than 50% 
covered) were forecast for areas of positive SD advection 
(negative relative vorticity advection). Note that clouds 
(greater than 50% covered) were assumed to be located only 
in areas of positive relative vorticity advection in this 
study. This is a slight deviation from Bittner’s forecast- 
ing procedures for clouds as noted below. 

7) To produce the second cloud cover forecast, steps 
1) through 4) are the same but step 5)is omitted since the 
SR prognosis was used as the advection field for the pur- 
pose of step 6) . 

8) To construct the third cloud cover forecast, steps 
l)through 4) are the same but step 5) is omitted since the 
unmodified Z prognosis was used as the advecting field 
(Bittner’s approach [1]) for the purpose of step 6). 

9) Steps 1) and 2) above were then applied to the 
next day’s ATS-1 phot ograph (current plus 24 hours) to 
produce a verifying nephanalysis . 

10) The verifying nephanalysis was separately put in 
register with the three forecasts. Each of the 351 grid 
points for each of the three forecasts was checked and 
compared with the corresponding grid point of the verify- 
ing nephanalysis. Grid points which were observed clear 
but forecast cloudy were denoted by a ”2”. Those which 
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were both forecast and observed cloudy were indicated with 
a Next, the points which were observed cloudy but 

forecast clear were denoted by a Finally the grid 

points which were not covered by clouds on the verifying 
nephanalysis nor forecast to be cloudy were represented by 
a "3". At this point the technique vrhich gave the best 
results could be readily determined by totaling the number 
of correct forecasts ("0” and **3") for each forecast day. 
Tables I, II and III give the raw scores for the advecting 
fields of the Z' PROG, SR PROG and Z PROG, respectively. 

The sixth column represents the number of correct forecasts 
while the seventh column denotes the number of incorrect 
forecasts . 

11) Finally, the National Meteorological Center’s 
threat score was computed for each of the forecasts. An 
explanation of the threat score is presented in Section VI. 

As indicated earlier, the manner in which the clouds 
were determined by this author is not strictly the same a? 
that used by Bittner [1]. In addition, he applied the rules 
of Appendix B in reverse to indicate some geometric detail 
in the cloud patterns, as spiral clouds, vortices, etc., 
whenever they were deemed necessary as determined by the 
SD’ prognosis. He also used various crossing-angle rules 
as a guide for forecast cloud tops, areas and amounts of 
precipitation, etc. Since this study was conducted to try 
and improve Bittner’s basic premise, these refinements were 
not included. 
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TABLE I. Verification of grid-point cloud forecasts using the Z’PROG as the 
advecting field. (Tabular values are number of grid points.) 
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TABLE II. Verification of grid-point cloud forecast 
advecting field. (Tabular values are number of grid 
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TABLE III. Verification of grid-point cloud forecasts using the Z PROG as the 
advecting field. (Tabular values are nunber of grid points.) 
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V. 



CASE STUDY OF MODIFICATIONS TO ANALYSES 



The synoptic time of 0000 GMT 10 April 1972 and its 
associated analyses and prognoses were randomly selected 
as an example of the procedures presented in the preced- 
ing section. Figures 1 through 18 graphically illustrate 
the following discussion. 

Figure 1 is the ATS-1 photograph of 2343 GMT 9 April 
1972 which was selected to be representative of the cloud 
cover for 0000 GMT 10 April 1972. The major cloud fea- 
tures (areas of ^ 50% cloud cover) are outlined for 
clarity . 

In general, SD patterns moving into the western sector 
of the verification grid were modified very little because 
the cloud patterns are not well-defined near 170®E in the 
ATS-1 pictures. This sometimes led to a large error in 
the verification of cloud forecasts west of 180®. There 
was usually good definition of the cloud features and, 
therefore, accurate modifications of the SD patterns in 
the extreme eastern section of the grid. These eastern 
SD features and associated cloud features sometimes moved 
out of the area of the verification grid during the time- 
span of the prognosis. Some of the centers of the SD 
analyses in this near North American coastal area were 
extensively modified. For those systems which did not 
appear within the verification area at the verifying time 
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of the prognosis, the overall goodness of the associated 
modifications and the resultant cloud forecasts are un- 
known. 

Figure 2 shows a nephanalysis (areas of ^ 50% cloud 
cover as derived from Figure 1, on a grid of the same 
scale . 

Figure 3 is the nephanalysis after it has been trans- 
ferred to a polar stereographic projection with the SD 
and SR fields superimposed. The SD pattern is illustrated 
with dashed lines and the SR field with solid lines. Note 
the strong SD low center at 43®N and 135'^W. Such a low 
value for the SD center suggests a cyclone at 500 mb. 

Since this center is embedded in an SR trough, as indicated 
by the solid lines, there is high probability of a closed 
center at 500 mb. However, the cloud patterns of Figure 1 
indicate only slight cyclonic circulation with convective 
clouds in the area. This suggests decreasing the positive 
vorticity advection (negative SD advection) considerably. 

Figure 4 shows the new or modified SD field. Note 
the drastic changes made to the SD zero line. In the lower 
left portion of the figure, the zero line was relocated to 
the poleward edge of the frontal cloud band. (Rule l.a 
Appendix B) . In general, SD troughs and ridges were re- 
drawn to coincide with major cloud boundaries. Note the 
considerable filling of the SD low center at 43®N, 135°W 
to decrease the negative SD advection. The center at 33°N 
180° was moved slightly northward to maintain the correct 
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perspective with the circulation center on the photograph. 
(See rule 2.e Appendix B.) 

Figure 5 is the modified SD analysis of Figure 4 
(dashed lines) superimposed on the original SD field of 
Figure 3 (solid lines). The amount and direction of move- 
ment of all centers and isolines were noted and used in 
the modification of the SD prognosis. 

Figure 6 represents the 36^hour SD prognosis and its 
modification (SD* PROG). As such, it represents the 
changes made to the SD analysis 24 hours earlier, as shown 
in Figure 5. As can be seen from the figure, this modifi- 
cation is somewhat subjective. 

Figure 7 shows the graphical addition of the modified 
SD prognosis (SD* PROG) to the SR prognosis in order to 
produce a new height prognosis (Z* PROG). This is the 
only additional step required over Bittner’s method. 

Figure 8 is the Z’ prognosis of Figure 7 superimposed 
on the SD* prognosis of Figure 6. Since the negative SD 
values correspond to positive values of vorticity, as 
shown in Appendix A, the areas of positive vorticity advec- 
tion (PVA) are readily apparent. These areas are carefully 
delineated and are represented in the figure as clouds. 

Note that major cloud boundaries are determined by the SD 
troughs and ridges as defined earlier relative to vorticity 
advection and as stressed by Bittner. As mentioned in the 
previous section, Bittner inserted cloud-pattern detail 
(i.e. cloud vortices, developing cyclones, etc.) at this 
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point. These modifications could be done by a knowledge- 
able synoptic meteorologist resulting in better continuity 
of cloud patterns and increased accuracy of cloud fore- 
casts. 

Figures 9 and 10 are the same as Figure 8 except that 
the vor t ici ty-advect ion areas are determined by using the 
SR unmodified Z prognoses as the advection fields, respec- 
tively, See procedures 7 and 8 of the previous section. 

Figure 11 is the verifying ATS-1 photograph with areas 
of ^ 50% cloud cover outlined. This picture also shows 
very little cyclonic motion in the northeast corner, justi- 
fying the decreased PVA of Figure 4, 

Figure 12 shows the boundaries of the ^ 50% cloud cover 
from the photograph of Figure 11. 

Figure 13 illustrates the verifying clouds as trans- 
ferred from Figure 12 and the PVA areas derived in Figure 
8 from the Z’ prognosis. As mentioned earlier, the date 
for this case study was chosen at random; it just happens 
to be a below-aver age date. See Tables I, II and III for 
comparisons of raw values for 0000 GMT 11 April 1972, 

Figure 14 represents the accuracy of the forecast. 

The 351 grid points of Figure 13 (every 2-1/2® latitude- 
longitude spacing) were examined and assigned to one of 
four categories. The grid points that were forecast to 
be cloudy but were observed to be clear (”2**) totaled 65, 
Those that were both forecast and observed cloudy (**0”) 
amounted to 141, Observed cloudy but forecast clear (”!**) 
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were 101 points. And, finally, those grid points which 

were both forecast and observed to be clear (**3”) totaled 

43. This amounted to 184 correct (**0 + 3”) and 166 (**1 + 

2**) incorrect forecasts, which yields a combined threat 
2 

score of .67. See discussion in Section VI and scores 
in Table IV. 

Figures 15 and 16 show the verifying clouds superim- 
posed on the PVA area derived from the SR prognosis and 
the scoring grid, respectively. The number of grid points 
with correct forecasts (”0 + 3”) amounted to 170 while the 
incorrect total (’*1 + 2**) was 181. The threat score for 
this method was .59. 

In general, the scores for all three methods in this 
case study were below their respective averages. See 
Table IV for the averages. However, the relative differ- 
ences between the three methods are of the same order as 
the overall averages of each. That is the SR method of 
determining the PVA areas yielded approximately the same 
score as Bittner’s technique (unmodified Z field) for 
determining PVA areas while the Z* (modified Z field) 
m ethod gave better than a 10% increase in the value of the 
combined threat score which was used to evaluate the cloud 
cover forecasts. 



Threat score will be explained in detail in the following 
section . 
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VI. VERIFICATION OF RESULTS 



The National Meteorological Center’s threat score sys- 
tem was used to measure the goodness of the forecasts. 

The score depends on the correct forecasting of cloudy 
(Tq ) and clear areas (T^) • They are calculated as follows 



0 = Those points both forecast and observed as 

cloudy . 

1 = Those points forecast clear but observed as 

cloudy . 

2 = Those points forecast cloudy but observed as 

clear . 

3 = Those points both forecast and observed as 

clear . 

Using this approach, a perfect forecast of clouds re- 
sults in a score of T^ = 1.0 and a perfect forecast of 
clear skies results in a score of T^ = 1.0. Combining 
these two perfect threat scores gives a total of 2.0 for 
an absolutely correct forecast over an area as follows: 



The results of the three forecasts prepared for each 
day are listed in Table IV. These scores were obtained 
using the above technique. A cloudy forecast score and 
a clear forecast score followed by a combined score are 





0 3+1+2 



3 



where : 



T ^ + « 1.0 + 1.0 = 2.0 

O' 0 
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presented for each day of the study period. The overall 
averages of each method are listed at the bottom of the 
table • 

The results of Bittner’s study of approximately the 
same area of the North Pacific Ocean, but for different 
dates, are presented in Table V for comparison. The un- 
modified column represents threat scores obtained without 
the use of satellite pictures to modify either the initial 
or prognostic SD fields. The modified column represents 
the forecasts obtained after the initial SD fields have 
been modified based on the satellite observed cloud pat- 
terns and the prognostic SD fields have been modified 
based on the modifications applied to the initial SD 
patterns. This latter column represents a 15% increase 
in cloud-cover forecasting accuracy over the unmodified 
column. Except for the additional refinement steps men- 
tioned earlier, this modified column is equivalent to the 
Z PROG column of Table IV. The additional refinements in 
forecasting cloud cover used by Bittner may account in 
large part for the relative accuracy differences between 
the two corresponding columns of Tables IV and V. How- 
ever, as noted earlier, these refinements may be easily 
added to the basic approach in determining the PVA-associa- 
ted cloud cover. 

The basic determination of the PVA areas by the three 
different methods used in this study are clearly comparable. 
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The results in Table IV indicate that the Z PROG and the 
SR PROG advection fields are about equal in their results. 
The Z’ PROG advection fields indicate about a 10% improve- 
ment over the approach suggested by Bittner [1]. Thus it 
may be argued that if Bittner’s approach is modified to 
the extent of including the Z’ (vice the Z) prognosis for 
the advecting field, before adding the required and obvious 
continuity of large-scale cloud patterns, the best 24-hour 
cloud-cover fields may be produced. 
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TABLE IV. A comparison of threat scores of Z'PROG, SR 
PROG and Z PROG cloud-cover forecasting techniques. 



Verification 
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.63 


. 55 


.16 


. 71 


4/7 


.76 


. 22 


.98 


.72 


. 18 


. 90 


.70 


.18 


.88 


4/8 


.67 


.29 


.96 


. 63 


. 23 


. 86 


.62 


.34 


.96 


4/9 


. 63 


.44 


1.07 


. 63 


.43 


1.06 


.59 


.38 


.97 


4/10 


.59 


.29 


.88 


.75 


. 21 


.96 


.47 


.18 


.65 


4/11 


.46 


. 21 


.67 


.43 


.16 


. 59 


.41 


.18 


.59 


4/12 


. 89 


. 15 


1.04 


. 53 


. 16 


.69 


.53 


.19 


. 72 


4/13 


.44 


. 32 


. 76 


. 43 


.34 


.77 


. 38 


. 32 


. 70 


4/14 


.27 


.24 


.51 


. 26 


.24 


.50 . 


.34 


.25 


.59 


Average 


.56 


.31 


.87 


.51 


.29 


. 80 


.49 


.29 


.78 
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TABLE V 



The results of Bittner's study [1] 



DATE 


UNMODIFIED 


MODIFIED 




1969 


0 


"o 


«= T 


■^0 


+ ^0 


= T 


10/3 


.34 


.66 


1.00 


.22 


.61 


. 83 


10/9 


.30 


.63 


.93 


. 33 


. 73 


1.06 


10/10 


.29 


.49 


.78 


. 34 


.66 


1.00 


10/13 


. 24 


.58 


.82 


.32 


. 66 


.98 


10/14 


. 14 


. 60 


.74 


.23 


.63 


.86 


10/15 


.28 


.59 


.87 


.32 


.,68 


1 . 00 


10/16 


.24 


.57 


.81 


.34 


.74 


1.08 


10/17 


. 35 


.68 


1.03 


.26 


.70 


.96 


10/18 


.24 


.57 


.81 


.36 


. 80 


1.16 


10/19 


.45 


.79 


1.24 


.40 


. 84 


1.24 


10/20 


.27 


.62 


.89 


.31 


.68 


.99 


10/21 


.32 


.65 


.97 


.34 


. 71 


1.05 


10/24 


.26 


.43 


.69 


.25 


.68 


.93 


10/25 


.38 


.61 


.99 


.27 


. 60 


. 87 


10/27 


.25 


.60 


.85 


. 43 


.78 


1.21 


10/30 


.26 


.35 


.61 


.40 


. 63 


1.03 


11/4 


.32 


.57 


.89 


.28 


.57 


.85 


Average 


.29 


.59 


.88 


.32 


. 69 


1.01 
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VII. CONCLUSIONS 



This study indicates there is merit in using the satel- 
lite-observed cloud patterns to modify the position and 
orientation of 500-mb relative vorticity patterns (the SB 
fields of FNWC) in order to improve 24-hour forecasts of 
cloud cover over oceanic areas. Bittner [1]> in an earlier 
study, came to the same conclusion. 

It appears that Bittner's basic approach can be improved 
by using both modified prognostic relative vorticity pat- 
terns and modified prognositic height fields at 500 mb. 
However, the refinements of cloud-pattern continuity de- 
scribed by Bittner should be included in the final 24-hour 
cloud pattern forecast for best results. 
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Figure 1. 2343 GMT 09 April 1972 ATS-1 satellite photo- 
graph (Subpoint: 2®S). Area of ^ 50% cloud cover bounded 
by O . 




Figure 2. Nephanalysis of ATS-1 satellite photograph 
(Subpoint: 2®S), 2343 GMT 09 April 1972. Area of ^ 50% 
cloud cover bounded by O . 
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Figure 3. FNWC 500~mb SD ( ) and SR ( ) analyses, 0000 

GMT 10 April 1972, superimposed on the ATS-1 nephanalysis 
of ^ 50% cloud cover (2343 GMT 09 April 1972). SD interval: 
30 m; zero SD line labeled. 




Figure 4. FNWC 500~mb SR ( ; ana moairiea bu ^ ) 

analyses, 0000 GMT 10 April 1972, superimposed on the 
ATS-l nephanalysis of ^ 50% cloud cover (2343 GMT 09 
April 1972). SD interval: 30 m; zero SD line labeled. 
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Figure 5. 500-mb FNWC SD ( ) and modified FNWC SD ( ) 

analyses, 0000 GMT 10 April 1972; SD interval: 30 m; zero 
SD line labeled. 




Figure 6. 36-hour 500-mb FNWC SD ( ) and modified FNWC 

SD ( ) prognosis, verifying at 0000 GMT 11 April 1972; 

SD interval; 30 m; zero SD line labeled. 
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Figure 7. Modified FNWC Z prognosis ( ) as derived from 

the FNWC SR (-•-•) and modified SD ( ) prognosis, verify- 
ing at 0000 GMT 11 April 1972: SD interval: 30 m; zero SD 

line labeled. 




Figure 8. Modified FNWC 500-mb Z ( ) and SD ( ) prog- 

noses with resulting areas of positive vorticity advection 
( ), verifying at 0000 GMT 11 April 1972 : SD interval: 

30 m; zero SD line labeled. 
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Figure 9. FNWC 500-mb SR ( ) and modified SD ( ) prog- 

noses with resulting areas of positive vorticity advection 
( C? ), verifying at 0000 GMT 11 April 1972: SD interval: 

30 m; zero SD line labeled. 




Figure 10. FNWC 500-mb Z ( ) and modified SD 

noses with resulting areas of positive vorticity 
( (p )» verifying at 0000 GMT 11 April 1972 : SD 
30 m; zero SD line labeled. 



( ) prog- 

advection 
interval : 
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Figure 11. 2358 GMT 10 April 1972 ATS-1 satellite photo- 

graph (Subpoint: 2°S). Area of ^ 50% cover bounded by 




Figure 12. Nephanalysis of ATS-1 satellite photograph 
(Subpoint: 2°S), 2358 GMT 10 April 1972. Area of ^ 50% 
cover bounded by • 
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Figure 13. Nephanalysis of ^ 50% cloud cover from ATS-1 
photograph (2358 10 April 1972) ( O ), and area of posi- 

tive vorticity advection, (((jv*.’.- ) derived from the modified 
36-hour FNWC 500-mb Z prognosis verifying at 0000 GMT 
11 April 1972. 




■ . . f • . 0 . • -y •'* \ • ■ . • ■ . • 



Figure 14. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb Z prognosis, verifying at 0000 
GMT 11 April 1972. Legend: 0-Clouds: forecast and observed- 

141 points; 1-Clouds: observed but not forecast-101 points; 
2-Clouds: forecast but not observed-66 points; 3-Clear: fore- 

cast and observed-43 points. 
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Figure 15. Nephanalysis of ^ 50% cloud cover from ATS-1 
photograph (2358 GMT 10 April 1972) ( C? ), and area of 

positive vorticity advection derived from 36-hour FNWC 
500-mb SR prognosis verifying at 0000 GMT 11 April 1972. 
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Figure 16. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb SR prognosis, verifying at 0000 
GMT 11 April 1972. Legend: 0-Clouds: forecast and observed- 

135 points; l^Clouds: observed but not forecast-105 points; 
2-Clouds: forecast but not observed-76 points; 3-Clear: fore- 

cast and observed-35 points. 
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Figure 17. Nephanalysis of ^ 50% cloud cover from ATS-1 
photograph (2358 GMT 10 April 1972) (O ), and area of 
positive vorticity advection derived from the 36-hour 
FNWC 500-mb Z prognosis verifying at 0000 GMT 11 April 1972. 
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Figure 18. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb Z prognosis, verifying at 0000 
GMT 11 April 1972. Legend: 0-Clouds: forecast and observed- 

127 points; 1-Clouds: observed but not forecast-110 points; 
2-Clouds: forecast but not observed-74 points; 3-Clear: fore- 

cast and observed-40 points. 
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Figure 19. The relationship between cloud patterns and the 
500-mb SD field, as presented by Bittner [1] and illustrated 
by Zeigler [8]: A. Surface wave vs. SD field. B. Occluded 
system with 500-mb flow and SD field. C. Mature cyclone vs. 
SD field. D. Dissipating cyclone vs. SD field. E. Polar jet 
axis and open-celled cumulus vs. zero SD contour. F. SD flow 
vs. cloud-associated PVA MAX. 
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APPENDIX A 



FNWC'S 500-MB SD FIELD AS AN APPROXIMATION TO 
THE RELATIVE VORTICITY FIELD 



The definitions and discussion below follow closely 
that in [ 1 ] . 

The 500-mb disturbance component is defined as 



aSD _ 

^^500 " ^ 500 ~ ^^500 “ ^ ^ 500 ^*^ 

o 



( 1 ) 



where : 



Z^Oo “ 500-mb height field 

= 500-mb residual (mostly wave numbers 1,2, 3, A, 

5 and 6) 

^^SOO “ 500-mb shortwave field (mostly wave numbers 
7 and above) 

a = degree-o f-smoo thing parameter 
6S = numerical grid length 

The right hand side of this equation is a Laplacian of the 
height field over the smoothing interval a . 

Geostrophic relative vorticity is defined as: 

C = gf. V (2) 

where: g = acceleration due to gravity 
f = Coriolis parameter 

By substituting Equation (2) into Equation (1), the follow- 
ing relationship is obtained: 

SD^qq = -6S fg“ f cda (3) 

o 
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Equation (3) shows that the SD field is about the same as 
the geostrophic relative vorticity field (with a sign re- 
versal) over the smoothing interval used in the analysis. 



Thus, for 
analogous 



all practical purposes, they can be considered 
for both analyses and prognoses. 
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APPENDIX B 



SUMMARY OF RULES RELATING CLOUD 
PATTERNS TO 500-KB SD ANALYSES 

Following is a summary of rules relating cloud patterns 
to the 500-rab analyses [1]. 

1) The SD zero line should be drawn: (See Figure 19 for 

a graphical illustration of some of these rules.) 

a) Along the poleward edge of the cloud band associated 
with the polar front. 

b) Along or just equatorward of the cloud shadow mark- 
ing the polar jet axis; along or just inside the 
poleward limit of transverse lines marking the sub- 
tropical-jet axis. 

c) Along the line between two major vorticies which 
are connected by a continuous frontal band. 

2) An SD low center should be located: 

a) Five to seven degrees upstream from an open frontal 
wave, or at the rear edge of the head of a cloud- 
identified positive vorticity advection maximum 
(PVA MAX) . 

b) Five degrees upstream from the center of a hook- 
shaped cloud pattern (rapidly developing cyclone). 

c) Within three degrees or less from the center of a 
complete spiral cloud pattern (mature cyclone). 

d) Directly over or slightly to the east of the center 
of an open vortex (dissipating cyclone). 
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e) 



Over the center of a strong vorticity center in 



the cold air. 

3) An SD trough should be drawn: 

a) Through the point where a frontal cloud band 
forks or fractures. 

b) Five to seven degrees upstream from a weak frontal 
wave (where a closed SD low is not justified). 

4) An SD ridge should be drawn: 

a) In an occlusion, to coincide with the downstream 
major cloud boundary or the change from C (cumuli- 
f orm/s trat iform) to MCO (s trat if orm/cirrif orra) . 

b) With anticyclonic curvature near the equatorward 
edge of a frontal cloud band. 

c) With anticyclonic curvature over the peak of a 
frontal wave pattern. 

5) The SD gradient should be maximized: 

a) Along the rear edge of a frontal cloud band. 

b) In the dry tongue of a hook-shaped vortex. 

c) At and poleward of the peak of a frontal wave. 

6) SD lines should be drawn: 

a) Parallel to lines of cumulus (CELLS OPEN) in the 
cold-air sector of a cyclonic disturbance. 

b) Parallel to cloud lines and bands in a dissipating 
vortex. 
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